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1. Introduction 

A vast literature on phenomenological descriptions of the rj—rf system has been published 
in the past ten years (T|. Yet, the 77 — 77' mixing angle alone is more than an effective 
parameter to be extracted from low energy data. Its peculiar value may indeed shed some 
light on the non-perturbative dynamics of the fundamental QCD theory and in particular 
on the axial U (1) anomaly. Needless to recall here why the subsequent parity (P) and time- 
reversal (T) violations constitute a major puzzle in the Standard Model for electroweak 
and strong interactions 0. 

To link this axial anomaly with the observed mass spectrum for the pseudoscalar meson 
nonet, alternative paths based on the chiral perturbation theory or the large number of 
colours limit have been proposed. Among them, the chiral perturbation theory at leading 
order in p 2 and 1/N C proves efficient once the typical 20% corrections expected from the 
flavour symmetry breaking are duly acknowledged. 

Within this rather simple framework, the 77(7/) masses are functions of the mixing angle 
9. In particular, the 77 — 77' mass ratio is not fixed by the theory but can only be optimized 
with respect to its experimental value for 9 ~ —27°. However, the corrections requisite to 
reproduce the measured value of this ratio raise the question of the systematic expansion 
to adopt. It appears that including the next to leading order in p 2 in the large N c limit is 
quite predictive and compatible with the data. Consequently, this approach requires the 
1/iVe-suppressed one-loop contributions to be small. In this Letter, we emphasize that the 
optimal value of the 77 — 77' mixing angle at leading order turns out to consistently damp 
out the quadratically divergent one-loop corrections to the 77 — 77' inverse propagator matrix 
and the 77' — » r/nir decay amplitude. 
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2. An effective theory at leading order in p 2 and 

If n quark flavours are massless, the fundamental Lagrangian of QCD displays a global 
U{n)i, <S> U(n)n invariance. In the large N c limit, N c being the number of colours, the 
effective Lagrangian which features this chiral symmetry at lowest order in p 2 reads || 

f2 r 



(2.1) 



where U is a n-by-n matrix transforming as U — > giJJg R . The mass matrix m for the 
light quarks transforms as U and its determinant is assumed to be real to ensure P and T 
invariance. In Eq.Q2.ip, the parameters with dimensions of mass scale respectively as 



/ oc N^ 2 , r oc AT,' 



(2.2) 



In the large N c limit, U(u)l <8> U(n)n has to be spontaneously broken into the maximal 
vectorial subgroup U(n)y if n > 3 0]. Consequently, U is a unitary field which can be 
expanded around its vacuum expectation value as a function of the Goldstone bosons nonet, 



U = 

In the case of three light flavours 
v V2K- 



l + iV2-- — + 

f f 



(2.3) 



3'l 



v / 2vr+ 
V2K° 



3 'I 



V2K+ 
V2K° 



(2.4) 



and the masses of the pseudoscalars can be easily extracted once m is diagonalized. Working 
from now in the isospin limit m u = = m, we obtain 



rm 



m lr = o 0" + ms 



(2.5) 
(2.6) 



and 



4m 2 K 



m 



ml 



-2\/2i 



K 




(2.7) 



with the octet-singlet flavour basis conventionally characterized by the following amount 
of strange/non-strange quarks in the meson wave function 



1 



r] ~ —= (uu + dd — 2ss) 
V6 

rf ~ —= (uu + dd + ss) . 
v3 

At this level, the masses of the physical pseudoscalar fields 




(2.8) 
(2.9) 

(2.10) 
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are only functions of the tt and K ones and vanish in the chiral limit m u = m^ = m s = 0. 
However, the measured mass of the rf around 1 GeV tells us that the axial U(l) has been 
broken by the dynamics of QCD itself ||. In the limit of a large number of colours within 
chiral perturbation, this explicit breaking is implemented through the one and only term 

§ 



f 2 ml 
8 4N r 



ln[/-lnC7 t 



1 



mlrf + 6 (tt 4 ) 



(2.11) 



which is l/-/V c -suppressed but p°-enhanced with regard to the effective Lagrangian ([2.1|). 
Accordingly, the 770 — t]q element mQ of the mass matrix ( |2.7D is corrected by the parameter 
mg so that the 77, r\ masses are not anymore fixed in terms of the it and K masses but are 
functions of the mixing angle 9, as displayed in Fig|l|: 



771.7 



77) 



rf 



4:171k 



Am 2 K 



ml + 2^2 [m\ 



tan 9 



ml - 2V2 (m 2 K - ml) cot I 



(2.12) 
(2.13) 



The resulting relation between physical quantities defined at lowest order 



tan 2 6> 



77)7 



(4m 



K 



ml 



(4m 



K 



mt 



mz 



11.4°) 



(2.14) 



is analogous to 



tan 9\y 



?))- 



77). 



w 



77), 



w 



mf 



(I' 



'w\ 



28.2°) 



In other words, the Gell-Mann-Okubo (GMO) mass relation m 2 



(2.15) 



8S 



I (4.m\ — ml) in the 

778 — 770 mass matrix ( |2.7D plays here the role of the isospin mass relation tti 2 ^ = m 2 v± 
in the W3 — Bq mass matrix of the Standard Model for electroweak interactions. The 
latter relation is known to be invariant under the unbroken custodial SO (3) of the Higgs 
potential; the former is invariant under the unbroken vectorial SU(2)f (g) U(l)y since the 
quark mass matrix m in Eq.([2.l|) transforms at most as a singlet and an octet of SU(3)v- 
A breaking of the GMO relation for mg 8 would require & (p 4 )0) terms like (mU^mU^ 
with m (g) m also transforming as a 27 under the vectorial flavour group. 



Surprisingly, even with the additional parameter mg, the masses of 77 and rj' cannot be 
fitted simultaneously ||. Indeed, taking away m K from Eqs.(2.12- 2.13| ), we easily obtain 



777; 



771: 



mt 



mi 



= tan (29th - 9) tan 9 (tan 29 th 
< tan 2 9 th = 2- Vs. 



-V2 



(2.16) 



In the safe m 2 — > limit, the resulting upper bound of 0.27 for the 77 — 77' square mass ratio 
is clearly at variance with the corresponding experimental value of about 0.33. 
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Figure 1: The r\ and rf masses as a function of their mixing angle from Eqs (2.12) and ( 2.13| ). We 
choose to work with £ [-f, +?] to avoid the renaming 77 — ► 77', 77' — > —77 at 9 = —j. If m Vy K 
are fixed at their experimental values, the measured 77 and 77' masses denoted by dots cannot be 
simultaneously reproduced at lowest order. 



Mass corrections of about 20%, as requested by Eq.( 2.16| ) to reproduce the observed 77—77' 



spectrum, drastically change the absolute value of the mixing angle derived in Eq.( 2.14| ). 
Indeed the physical mass of the 77 and the octet mass m^g turn out to be numerically close, 
within a few percent. Therefore, any departure of lowest order 77 mass from its physical 
value is enough to produce a major modification of the angle 6 extracted with the help of 
Eq.( [2.12 ), as illustrated in Fig.Q. So, a determination of the mixing angle at lowest order 



is sensible only if its value is stable with regard to 1/N C and chiral corrections. In this 
respect, any enlarged symmetry beyond the custodial one is welcome to tame the quantum 
corrections. For example, a parity-conserving local SU(3)l ® SU(3)r extension of the 
SU(2)l (8) U(l)y electroweak gauge symmetry H covers the custodial SO(3) and would 
imply 

tan#vK = — 7= (6 w = -30°) (2.17) 

in pretty good agreement with the on-shell absolute value of the weak mixing angle already 
introduced in Eq.( 2.15| ). 



In Eq.( |2.1| ), the canonical kinetic term for the tt field has a global SO(9) invariance. 
Both the vectorial 5C/(3)-breaking in Eq.( |2.l[) and the axial [/(l)-breaking in Eq.( 2.1l| ) 



already violate this symmetry at the level of the terms quadratic in the meson fields. Yet, 
for particular values of the angle 6, remnants of SO(9) may survive at this level; they 
correspond to the two mass degeneracies displayed with dashes in Fig.|l]: 
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if e = e id with 

tan ^=_L (e id = +35.3°), (2.18) 

the physical rf ~ (uu + dct) is degenerate in mass with the pions Q while rj ~ — ss. 
Note that the negative value 0^ = —54.7° corresponding to the other convention with 
the ss component singled out, namely r] ~ (uu + and rf ~ +ss, is outside 
the interval [-f , +§] (see Figj|). The ideal mixing obtained from Eq. (|2.7] ) , i.e., for 
uiq = 0, is relevant for the vector meson mass spectrum on which the axial U(l) 
anomaly has no effect, but totally unrealistic for the pseudoscalar one. 

If = 9 p h with 

tan 9 ph = 0= (9 ph = -19.5°) , (2.19) 

the physical r\ ~ (uu + dd — ss) is degenerate in mass with the kaons while 
rf ~ ^= (uu + dd + 2ss). Here, this sensible value for the mixing angle is called 
phenomenological since it has been extensively used to study hadronic B decays and, 
in particular, to explain the striking suppression of B — > Kr\ with respect to B — » Krf 



H if penguin diagrams dominate these processes [10]. It is also quite popular be- 
cause the associated quark components are easy to remember and to handle in a 
phenomenological quark-diagram description of the decay amplitudes according to 
their SU(3) properties. 



We have no simple mass degeneracy for the case of 6th already introduced in Eq.( 2.1€| ) 
but note that the three angles of peculiar interest are related through 

tan 26 th = tan (9 ph - 6 ld ) (6 th = -27.4°) (2.20) 

with, quite incidentally, 0th ~ @w if the weak mixing angle turns out to be negative as 
predicted by some unification theory. 

With respect to possible enlarged symmetries covering the custodial SU(2)i®U(1)y, we 
observe that the mass degeneracies mw = and m v = mx correspond to the breaking 
patterns £0(9) -> SO (A) <g> SO (A) and SO (9) -> SO (3) SO (5), respectively. These 
patterns for Q id and 9 p h can be understood from the fact that 50(9) group admits 517(2)® 



SU{2) ® 5p(4) or, equivalently, 50(4) ® 50(5) as a maximal subgroup [11]. However, 
such enlarged symmetries are explicitly broken at the level of the full effective theory and 
thus accidental. Consequently, the finite value of the Oid and Q p h mixing angles should 
not be protected against (quadratically) divergent quantum corrections. The fact that the 
relations pl8| ) and (pT9| ) are not natural can easily be confirmed through the following 



one-loop computation. 

3. One-loop corrections to the r\ — inverse propagator matrix 



The unification value ( |2.17 ) for the observable weak mixing angle 9\y can most easily be 



derived by requiring the one-loop fermionic contribution to the Z — 7 mixing diagram to be 
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finite [12]. In the same spirit, let us impose the cancellation of the quadratically divergent 
one- loop corrections to the rj — rf mixing angle 6. 



In order to compute these corrections, we need now to expand U up to the order 7r 4 , 

U = l + J2ak[iV2j) . (3.1) 



The parameter a\ may be absorbed into the definition of / while the even coefficients are 
fixed by the unitarity condition [13] 

oi = 1, a 2 = -, a 3 = b, a 4 = 6--, ... (3.2) 

l o 
with 6 an arbitrary parameter. For b = g, we recover the standard form 

U = exp(i^\ (3.3) 



also suited for an octet of pseudoscalars p4 |. But as shown in ref.[|lj|, any other value 
of b gives rise to the same T matrix when all external lines are put on the mass shell. 
Yet, one-loop corrections from the kinetic part of the Lagrangian ( [2.1| ) induce in principle 
a momentum-dependent rj — rj' mixing term which thus has to be taken off-shell. Again 
by analogy with the scale dependent Z° — 7 mixing induced at one-loop in the Standard 
Model, let us therefore introduce the propagator formalism [16]. 

If we denote by —iA XlX2 (p 2 ) with Xii X2 = f]i v' the one-loop contributions to the cor- 
responding two point functions, the inverse propagator matrix £ can be parametrized as 
follows 

S w = (1 + Z v ) (p 2 - m 2 ) + 5m 2 - A m (p 2 ) 
Z v , v , = (1 + Z v ,) (p 2 - m 2 ,) + 5m 2 , - (p 2 ) (3.4) 
S w / = bm 2 m , - A m i (p 2 ) . 

The last relation in Eq.(3.4) takes into account the fact that r] and rj are decoupled at 
tree-level, but leaves open the possibility for the one-loop induced mixing to depend on p 2 . 
Imposing the normalization of the kinetic part of to be canonical and the physical 

masses m Xi to be the poles of the propagators, we identify 

and 

5m li = A XiXi( m2 Xi ) ( 3 - 6 ) 
where the prime denotes the derivative with respect to p 2 . From a one-loop computation, 
we obtain the following quadratic dependences on the ultraviolet momentum cut-off A: 

A 2 

Z v = 3 [(3 - 206) + (46 - 1) cos 26} ? 

(4tt/) 

A 2 

Z v , = 3 [(3 - 206) - (46 - 1) cos 26} ? (3.7) 

(4tt/) 
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and 



6 (m 2 + m 2 ,) 
5 (m 2 m 2 t) 



-2 (2m 2 K + ml] 



A 2 



{MY 



-6m 2 (2m|- 



A 2 



{WY 



(3.8) 



with 



A 



rftf 



[3(46 
+4\/2(26 



1) 



m 2 ) cos 20 I ■ 



(3.9) 



(/r) ^ :■!( i// - l)p 2 + 2(1 - 8b)m 2 K + 2(26 - l)m 2 ] sin20 

A 2 
J (4vr/) 2 ' 

Here, the pseudoscalar masses mx,n and the mixing angle 8 are parameters associated with 
the lowest order Lagrangian defined by Eqs (2.1) and ( 2.11] ). In particular, m 2 , has been 
taken away with the help of the relation 

ml = ! (l - 2\/2cot26>) (m 2 K - m 2 ) . (3.10) 

In general, the one-loop quadratic divergences can be absorbed by a redefinition of the 
parameters in the & (p 2 ) Lagrangian. Indeed, the corrections quadratic in the cut-off can 
be identified with the d = 2 pole in dimensional regularization. Here, a full cancellation of 
the G (p 2 , divergent correction ([O]) to the mixing requires 



tan 28 (p 2 ) 



4^2(26-1) 



m 



K 



(3.11) 



3(1- 46) p 2 + 2 (86 - 1) m\ + 2 (1 - 26) m 2 ' 

Depending on the parameter b, the mixing angle defined in Eq.( |3.1l] ) is not a physical 
quantity. The only way to get rid of the b-dependence is to choose p 2 = 2m\. At such a 
momentum consistently located between the r\ and rf masses, Eq.( |3.11 ) then provides us 
with an effective mixing angle 8 defined at the QCD scale m 2 , : 



tan 28 (ml) 



-2V2 (m 



mZ 



-25.8 C 



(3.12) 



(2m\ + ml) 

We note that the same expression for an on-shell mixing angle 9 can be obtained by simply 
fixing 6 = \ to cancel the momentum dependence in Eq. (]3.9| ) . This value of the parameter 
b, which suggests the other significant form 



U 



1 + vff 



(3.13) 



only suited for a full nonet of pseudoscalars [ 13 1 , ensures ^-independent wave- function renor- 
malizations, i.e., = Z^t in Eq.(3.7). As a consequence, the only chiral invariant mass 
operator that would absorb any divergent rjg — ??o rotation at G (p 2 , l/N c ) is proportional 
to 

f 



! —r (mU ] - Um^ (in U - In = (2m 2 K + m 2 ) jjg - 2^2 ( 



m\ 



VOVS 



+ t 



IT 



(3.14) 
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in full agreement with Eq. fl3.8j ) and Eq.( 3.12| ). So, the parity-conserving global SU(3)i (g> 



SU(3)r plays here the role of the enlarged symmetry which covers the custodial SU(2)j <S> 
U{l)y- Indeed, Eq,( |3.14j ) tells us that the chiral symmetry of the full effective theory 
selects in a natural way one negative value (9) for the Tj — rf mixing angle, without spoiling 
the GMO mass relation for mj 8 . 

As already anticipated from the explicit breaking of the accidental symmetries SO (4) <S> 
SO (A) or SO (3) (8> SO (5) at the level of terms quartic in the meson fields, neither 9^ nor 
9 p h are protected against A 2 quantum corrections. On the contrary, Eq,( 3.12j ) tells us that 



the angle 9th which optimizes the 77 — 77' mass ratio at lowest order might be natural in 
the safe limit m 2 — ► 0. In the fundamental theory (i.e., QCD), the corresponding limit 
777"u,d ~~ * would, in principle, solve the so-called strong CP problem. This rather intriguing 
link evidently calls for further investigations. 



4. One-loop corrections to the 77' — > r/nn decay amplitude 

For the purpose of computing a 6-independent one-loop correction involving the 77 — 77' 
mixing, let us now consider a physical process with on-shell 77 and 77' states. 

4.1 Tree-level amplitude 

The tree-level amplitude for the 77' — ► rj-rnr decay reads 

A (77' -» 777rvr) = J2 2 ( 2 ^2 cos 29 - sin 2fl) Q - bj (m 2 + m 2 , + 2m 2 ) 
+8 (2V2 cos 29 — sin 29) lb ) rrh 



+4V2 (cos 2(9 - v / 2sin26») (b - ^\ ml 



(4.1) 



where m^, m„/ and stand now for the physical masses since they come from the mo- 



mentum dependence induced by the kinetic term in (|2.l| ). In Eq.(4.1), the second term 
proportional to r is due to the mass term in Eq. (|2.1| ) and the third one arises from the 
anomalous part given in Eq. ( |2.1lD . With the help of Eq. (ETlC|) , we eventually recover the 
well-known result that the tree-level amplitude 



2 

A (7/ -» ^vrvr) = ^ ^2 V2 cos 29 - sin 2#) 



3/ 

vanishes if 9 = 9id and is by far too small to reproduce the measured decay width. 



(4.2) 



4.2 One-loop amplitude 

The one-loop corrections to the process 77' — * rjmr are associated with the diagrams given 
in Fig.pl. 
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Figure 2: One-loop topologies for the r)' — > 777T7T decay amplitude. 



The first topology corresponds to the corrections of the inverse propagator given in 
section |3|. The second one involves 7r 6 vertices and thus requires the introduction of the 
next two coefficients in the development (fO) , namely 



o 5 = c 

b 2 b 
ae = C + J-2 

As a result, the A 2 -correction to the decay amplitude is given by 

,2 



1 

16' 



(4.3) 



5 A (7/ 



777 ' 

V nir) = 4-£ cos d 20 

r 



( tan 261 + v^l f tan 2 20 + — ^ tan 20 + - 
V / V 4a/2 2 



3 

+- 



m 



( tan 20 + 



4 mjj- — m 2 \ 



2V2 



tan 2 20 



4^/2 
A 2 
(4tt/) 5 



(4.4) 



This correction is independent of b and c, as it should for any physical quantity, and can 



be reproduced using the output of FeynRules [17] and Feynarts p8| . 

If we consider again the limit m 2 <C rn 2 K , we conclude from Eq. Q4.4D that the optimal 
value 6th given in Eq.( 2.20| ) for the r\ — rf mixing angle indeed damps out the quadratic 
dependence on the ultra-violet momentum cut-off A, as anticipated from Eq.(3.12). 



5. Comments and conclusion 

In the past, alternative ways to merge the large number of colours limit into the chiral 
perturbation theory have been used to study the rj — rj' system. In particular, the combined 
expansion 



V 



'(*) 



(5.1) 



has been advocated in ref.[19|. In this Letter, inspired by the pseudoscalar mass spectrum, 



we rather follow the approach of ref.[20| where the leading term in the 1/N C expansion is 
retained at each order in p 2 . At the effective level, this implies the hierarchy 



namely 



0(p°,l/N c ) >0{p 2 ,O) >ff(p 4 ,0) 
0(p 2 ,l/N c ) « €?(p 4 ,0) 



(5.2) 
(5.3) 
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with the large N c limit denoted by a zero as in Eqs.(2.1) and fl2.ll ). It amounts to remove 
the double trace term (3.14) as well as {d^UU^ (d^U^U^ in the Lagrangian, and to neglect 
the quadratic one-loop divergences which would renormalize them. The rf — > rnrir decay 
amplitude and the 77 — 7/ mass ratio are known to require sizeable corrections beyond the 
G (p 2 , 0) approximation and can thus distinguish between the two working hypothesis ( [5.1| ) 
and (1^). In ref.Q and ref.g, the G (p 2 , 1/N C ) contributions were invoked for the decay 
amplitude and the mass ratio, respectively. On the contrary, in ref.[22] and ref.pO|] the 
G (p 4 ,0) contributions were favoured for these physical quantities, respectively. 

At G (p 4 ,0), the full set of corrections allows us to naturally reproduce the observed 
rj — rj' mass spectrum. They do not fix by themselves the value of the mixing angle 9 but 
imply a splitting among the pseudoscalar decay constants [20]. In particular, the measured 
SU (3)-splitting between tt and K decay constants, 

f -f = l + e (5.4) 

Jit 

with e = 0.22 ± 0.01 of the order of {m 2 K — m 2 ) /lGeV 2 , provides a rather interesting link 
between our present work on the rj — rj' mixing and the so-called two-mixing-angle scheme 
high-lighted in ref.j||. Indeed, the equations 

0-0 2V ^f 

9 = 9+^e (5.5) 

relate the universal mixing angle which diagonalizes the octet-singlet mass matrix (after 
renormalizing the meson fields) to the #8,0 angles associated with the octet-singlet decay 
constants 

A = (i + |)/* 

/o = (l - |) fx- (5.6) 



3 

At G (p 2 ,0), e = and 9% = 60 but 9 cannot be determined. Yet, in this Letter, we have 
explicitly checked that the mixing angle 



-tan" 1 v^- -27° (5.7) 



vth-- 2 

which optimizes the rj — rj' mass spectrum at lowest order is protected against quadratic 
one- loop divergences in the safe m 2 — > limit. This result vindicates the approach based 
on Eq.( |5.2| ) since 9th is quite consistent with the physical mixing angle 

0«-(22±l)° (5.8) 

directly extracted from the anomalous J/^ — > rj(rj')j decays pCf. Indeed, higher order 



corrections are typically of the order of 20%, as nicely illustrated in Eq.(5.4). In conse- 
quence, 9g ~ —34° and #0 ra —10° within our specific momentum expansion supplemented 
by a large N c limit. However, any physical process only evaluated at the lowest order in 
the chiral expansion should rely on Eq. (|5.7| ) if it involves on-shell or off-shell r\ (7/), as it is 
the case in r\ (7/) — > 77 or in Kl — ► {rj, rj 1 ) — ► 77 decays, respectively. 
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